Correlation amplitude and entanglement entropy in random spin chains 
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Using strong-disorder renormalization group, numerical exact diagonalization, and quantum Monte Carlo 
methods, we revisit the random antiferromagnetic XXZ spin- 1/2 chain focusing on the long-length and ground- 
state behavior of the average time-independent spin-spin correlation function C(/) = In addition to the 
well-known universal (disorder-independent) power-law exponent T] = 2, we find interesting universal features 
displayed by the prefactor v = Do/3, if / is odd, and v = X>e/3, otherwise. Although Vo and De are nonuniversal 
(disorder dependent) and distinct in magnitude, the combination Do + De = — 1/4 is universal if C is computed 
along the symmetric (longitudinal) axis. The origin of the nonuniversalities of the prefactors is discussed in the 
renormalization-group framework where a solvable toy model is considered. Moreover, we relate the average 
correlation function with the average entanglement entropy, whose amplitude has been recently shown to be 
universal. The nonuniversalities of the prefactors are shown to contribute only to surface terms of the entropy. 
Finally, we discuss the experimental relevance of our results by computing the structure factor whose scaling 
properties, interestingly, depend on the correlation prefactors. 
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I. INTRODUCTION 
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Random low-dimensional quantum spin systems have been 
intensively investigated recently. The interplay between dis- 
order, quantum fluctuations, and correlations generates low- 
temperature phase diagrams with exotic phases.' In this con- 
text, one of the most investigated systems is the random an- 
tiferromagnetic (AF) quantum XXZ spin- 1/2 chain, whose 
Hamiltonian reads 



cr(0 = c?^(0 = (-i)' 



+ 1 



+ 1) 



(1.1) 



in which / labels the chain sites, S, are the usual spin- 1/2 op- 
erators, 7, 's are positive uncorrelated random variables drawn 
from a probability distribution Pq (J), and A,'s are anisotropy 
parameters, also random uncorrelated variables. 

The clean system, Jj = 1 and A,- = A, is a Tomonaga- 
Luttinger liquid for — 1 < A < 1 , with well-known asymptotic 
ground-state correlation functions, ^'^ 

Cf (/) = (5f5-;V/) = (-I)'f/-^'- -F/-^^-'/^^- , (1.2) 
C? (0 = {StSU^ = (-l)'A/-'/^c _ _L_ ^ (1.3) 

as I oo. The clean-system exponent is^ r|c = 1 — 
(arccos A) / 7t. At the "free-fermion" point A = 0, the pref- 
actors of the leading terms are known exactly,^'^ being given 
by A = 1/ (27i2) and F w 0.14709.'' Away from this point 
(|A| < 1), analytical forms for A and F were derived by 
Lukyanov and Zamolodchikov^'^ and checked numerically 
later on.'^ Furthermore, the constant F, evaluated numeri- 
cally in Ref. 9, is at least 1 order of magnitude smaller than 
F . At the isotropic point A = 1, irrelevant operators become 



(1.4) 



For A > 1, a spin gap opens and the system enters an anti- 
ferromagnetic Ising phase; otherwise, for A < — 1, the chain 
becomes a gapped Ising ferromagnet. 

Disorder strongly modifies the behavior in the clean criti- 
cal regime. It was shown that even the least amount of dis- 
order in Ji destabilizes the Tomonaga-Luttinger phase, pro- 
vided - 1 /2 < A,- < 1 . For - 1 < A; < - 1 /2, a finite amount 
of disorder is required to destabilize the clean phase. The 
low-energy behavior of the random AF spin- 1/2 chain then 
corresponds to a random-singlet phase, characterized by ac- 
tivated dynamical scaling with a universal "tunneling" expo- 
nent x|/ = 1/2, i.e., length (/) and energy (il) scales are re- 
lated through i2 ^ exp(— Z*''), irrespective of Pq (7).'^ More- 
over, the transverse and the longitudinal mean spin-spin cor- 
relation functions decay as a power law ~ X)l^^ for large dis- 
tances, both with the same universal exponent r| = 2.'^"'^ The 
mean value of the correlation function is dominated by rare 
widely separated spin pairs coupled in strongly-correlated sin- 
glet states. The remarkable fact that all correlations {xx, yy 
and zz) decay with the same exponent, irrespective of A, can 
be ascribed to the isotropy of the singlet state. In contrast, the 
typical value of the correlation function decays as a stretched 
exponential ^ exp(— /^). These results were obtained by us- 
ing the most successful theoretical tool to investigate such 
systems, the real-space strong-disorder renormalization-group 
(SDRG) method, first introduced in Refs. 14 and 15. 

The main idea behind the SDRG method is to gradually 
lower the energy scale by successively coupling the most 
strongly interacting spin pairs into singlet states. At each 
step of the renormalization transformation, one such pair is 
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decimated out of the chain, and its remaining neighboring 
spins become connected by a weaker renormaHzed coupling 
constant, calculated within perturbation theory. Thus, in this 
framework, the ground state can be viewed as a collection 
of "noninteracting" singlets formed by arbitrarily distant spin 
pairs. Although this description is not strictly exact, spin pairs 
do couple in states arbitrarily close to singlets.'^ 

Recently, efforts to compute certain numerical prefactors 
on disordered systems have been made. Fisher and Young 
have shown that the end-to-end correlation amplitude of the 
random transverse-field Ising chain at criticality is nonuniver- 
sal because of some high-energy small-scale features that are 
not treated correctly by the SDRG method. It is reasonable to 
expect that the same holds for bulk correlations. Indeed, no 
sign of universality was found in the correlation amphtude of 
the random XXZ chain. ^^-''^'^o 

Refael and Moore,^' on the other hand, have considered 
the mean entanglement entropy 5'(/) (for a recent review, see 
Ref. 22) between two complementary subsystems A (of size I) 
and B. Similarly to the clean system,^^'^^'^^ they have shown 
that S{1) = b + (y/3) In /, diverging logarithmically with the 
subsystem size. More interestingly, the prefactor y is uni- 
versal for a large class of systems governed by an infinite- 
randomness fixed point, namely, the random transverse-field 
Ising chain at the critical point and the spin- 1/2 random an- 
tiferromagnetic XXZ chain. Later, this amplitude was shown 
to be universal for a broader class of systems governed by 
infinite-randomness fixed points: the random ^-state Potts 
chain and the clock chain^^ and the random antiferromag- 
netic spin-S' chain at the random-singlet phase. ^^'^^ Moreover, 
it has also been shown that this amplitude is also universal in 
a large class of aperiodic chains. In the renormalization- 
group sense, and following Fisher and Young, Refael and 
Moore^' argued that the nonuniversalities of the correlation 
amplitudes are related to inaccuracies of order of the lattice 
spacing in the location of the effective spins. Such errors can 
only contribute a surface term b to the entanglement entropy, 
and therefore, its prefactor y should remain universal. No- 
tably, this should explain why all those other models display- 
ing a random-singlet-hke ground state show universal entropy 
prefactors. 

Our aim in this work is to further explore the issue of uni- 
versality in the behavior of ground-state correlation functions 
in random antiferromagnetic XXZ chains and make some di- 
rect links between spin correlations, structure factor, and en- 
tanglement entropy. We first calculate exactly, within the 
SDRG framework, the numerical prefactor v of the mean 
correlation function (S, • S,+/) in the limit I ^ °°, by relat- 
ing it to the distribution of singlet-pair bond lengths in the 
ground state. Surprisingly, it turns out to be universal and 
equals = — 1/4, if / is odd, and Ve = 0, otherwise, be- 
cause the "noninteracting" singlets can only be formed be- 
tween spins separated by an odd number of lattice sites. Nat- 
urally, this result is an artifact of the (perturbative) SDRG 
scheme, as shown by exact diagonalization (ED) studies of 
the XX model '^''^ (in which A, = 0, V;) and quantum Monte 
Carlo (QMC) calculations applied to the isotropic Heisenberg 
model^" (in which A, = 1, Vi). Nevertheless, as we show 



from ED calculations, in the XX limit, the long-distance be- 
havior of the longitudinal mean correlation function (5';5'j^/) 
is shown to be very well described by this renormalization- 
group prediction, while the transverse mean correlation func- 
tion {SJS'j^i') exhibits two distinct prefactors, and for 
odd and even /, respectively [see Eq. (3.3)]. 

Furthermore, we explicitly show that the mean entangle- 
ment entropy is directly related to the bond-length distribution 
of the singlet pairs and, therefore, directly related to the corre- 
lation function. This is interesting because it links a pairwise 
quantity (correlation) with a blockwise one (entropy). 

Since such a relation arises in the scenario of "noninteract- 
ing" spin singlets, in which nonuniversal effects are omitted, 
we introduce a toy model in which the correlations between 
different singlet pairs can be treated exactly. From the toy 
model, we gain some insight into the microscopic nature of 
the random-singlet phase and quantify the role of the interac- 
tions between the singlets. 

The picture emerging from our analytical results, and con- 
firmed by our own ED and QMC calculations, is the follow- 
ing. At long-length scales, the chain can be recast as a col- 
lection of noninteracting effective-spin singlets sharing strong 
pairwise correlations. These effective spins are clusters of 
original spin variables whose number depends on the details 
of the coupling constant distribution. With respect to the origi- 
nal spin variables, the singletlike correlations "spread" among 
the spins in the cluster, an effect which leads to the nonuniver- 
sal high-energy contributions discussed in the literature. 
Here, we go further by quantifying these contributions within 
the exactly solvable toy model. Interestingly, whenever the 
correlation is computed along a symmetry (z) axis, it equals 
the corresponding -1/4 singlet contribution, i.e., summing the 
correlations between all pairs of spins sitting at different clus- 
ters gives -1/4. For the mean correlation function, the result is 
that the combination of prefactors +\)^ = — 1/4 is univer- 
sal. For correlations along a nonsymmetry (x or y) axis, not 
only the prefactors are nonuniversal but also the combination 

+ ^e- This points to the importance of symmetry for the ob- 
served universality, a feature which was absent from the previ- 
ously considered models.'^ Finally, effective spins contribute 
only nonuniversal surface terms to the entanglement entropy, 
as expected. When one traces completely one of the spin clus- 
ters, its contribution to the entanglement entropy is the same 
as that of a singlet pair; only when the boundary between the 
subsystems is crossed by one of the clusters does that cluster 
contribute a nonuniversal term. Therefore, the entanglement- 
entropy prefactor is universal regardless of the existence of a 
symmetry axis. 

As a supplement, we compute both analytically and numer- 
ically the static structure factor S{q), which can be probed 
by neutron scattering experiments. Interestingly, we show 
that S (q) is dictated by disorder in the small-^ limit, namely, 
S{q 1) = k|^|. This is a consequence of two facts: (i) the 
decay exponent r| = 2 being universal, and (ii) the magnitude 
of the numerical prefactors Vo and Dg being different. More- 
over, we find that K = —n^{Vo +'Ue)/3, which implies that 
K = 7t^/12 is universal for5 computed along a symmetry axis. 
On the other hand, the behavior near the AF peak ^ = 7t is 
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dominated by the characteristic divergence of the clean sys- 
tem. However, the true divergence at ^ 7t is suppressed by 
disorder and the peak width is broadened. Since there is no di- 
vergence in the case of S (q) along the z axis in the XX model, 
disorder universally determines its behavior near the AF peak, 
i.e., s~ = 71 - e) = 7t -K|e| for e ^ 1. Only for q « n/2 is 
the clean-system behavior {q) = \q\ found. 

The remainder of this paper is as follows. We derive the 
universal SDRG expression for the mean correlation function 
in Sec. II, reporting our numerical analyses in Sec. III. Section 
IV discusses an exactly solvable model that yields instructive 
results on the origin of the universal behavior of correlation 
functions. In Sec. V, we derive the entanglement entropy and 
relate it to the distribution of singlet lengths and to the corre- 
lation function. We discuss the experimental relevance of our 
results by computing the structure factor in Sec. VI. Finally, 
we make some concluding remarks in Sec. VII. 



II. MEAN CORRELATION FUNCTION IN THE 
STRONG-DISORDER RENORMALIZATION-GROUP 
FRAMEWORK 

We start this section with a brief review of the SDRG 
method, followed by the derivation of the mean correlation 
function. 



A. Strong-disorder renormalization-group method: A brief 
review 

The main idea behind the SDRG method is to reduce the 
energy scale by integrating out the strongest couplings and 
renormalizing the remaining ones. In the present case, one 
locates the strongest coupling constant Q. = max {7,}, say, 
J2, and then exactly treats the two-spin Hamiltonian Hq = 
D. [S^S^ + S'-jS'^ + A25'|5'3) , considering Hi = // — //q as a per- 
turbation.'^ At low energies, spins S2 and ^3 "freeze" into a 
(nonmagnetic) singlet state, with the result that they can be 
effectively removed from the chain, provided that the neigh- 
boring spins and 54 are now connected by a renormahzed 
coupUng constant 



7 = 



(i + A2)n 



(2.1) 



calculated within second-order perturbation theory. The 
anisotropy parameter is also renormalized to A = AiA3(l + 
A2)/2. Note that / is smaller than either Ji, Jj, or £2, leading 
to an overall decrease in the energy scale. After the decima- 
tion procedure, the distance between and S4, which are now 
nearest neighbors, is renormalized to 



l = h+l2 + h 



(2.2) 



with /, defined as the distance between the spin 5, and its near- 
est neighbor to the right. The SDRG decimation scheme is 
illustrated in Fig. 1. 

Clearly, as the SDRG scheme is iterated and the energy 
scale Q. is reduced, the distribution of effective coupling 




J, A 



Figure 1: (Color online) Schematic decimation procedure (see text). 



constants Pj{J;Q.) is modified. Searching for fixed-point 
coupling-constant distributions {J;Q.), Fisher found that 
there is only one regular stable fixed point 



l-a 



(2.3) 



in which Q{x) is the Heaviside step function and with a = 
— l/ln£2 (we set the initial energy scale £2o to 1). This 
has been named an infinite-randomness fixed point (IRFP) 
since, as £1 0, the distribution becomes infinitely broad, i.e., 
VVar J/ J 00, where Var J and / are the variance and mean 
value of the coupling constants, respectively. Thus, the pertur- 
bative decimation procedure becomes more and more precise 
along the flow because the probability that both Ji and ^3 are 
much smaller than J2 increases as the energy scale is lowered. 

A useful quantity to be calculated is the fraction «n of "ac- 
tive" (not yet decimated) spins at the energy scale £2. It is 
obtained from the rate equation 



d«n = 2nnPj (J = £2; £2) d£2 



(2.4) 



where 2naPj {J = £2; £2) d£2 is the fraction of decimated spins 
when the energy scale is lowered from £2 to £2 — d£2. Hence, 
close to the fixed point. 



1 



(2.5) 



Equation (2.5) directly gives the low-temperature magnetic 
susceptibility x{T). One iterates the SDRG procedure until 
the low-temperature scale T . Spin pairs decimated at high en- 
ergy scales £2 S> r are "frozen" into singlet states, and thus 
their contribution to the magnetic susceptibiUty can be ne- 
glected. As the fixed-point distribution is very broad, all cou- 
plings between active spins are very weak compared to T , and 
the active spins can be considered as essentially free spins, 
each contributing a Curie term to the susceptibility.'^ There- 
fore, 



1' 



TlJ 
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1 



rin-'r 



(2.6) 



The low-energy modes are also given by Eq. (2.5). These 
modes are polarizations of widely separated weakly coupled 
singlet pairs, decimated at the energy scale £2 for which the 
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mean distance between spins was I n^^. Thus, the energy 
cost £2 to break a singlet of length / is 



£2-exp(-Z'f) 



(2.7) 



in which \|/ = 1/2.'^ This unusual exponential relation be- 
tween £2 and / is named "activated" dynamical scaling and 
\|/ has been dubbed the "tunneling exponent." 

The scaling behavior of the mean correlation function C (/) 
is cleverly obtained when one realizes that typically two very 
distant spins are not in a singlet state and thus are only weakly 
correlated. On the other hand, some rare and arbitrarily sep- 
arated spin pairs that were decimated together are strongly 
correlated and hence dominate the long-distance behavior of 
C{1). Therefore, the mean correlation function must be pro- 
portional to the total number of spin singlets decimated at the 
length scale /. Since the probability of decimating a spin pair 
is proportional to the probability that both spins have not been 
decimated yet, it follows that 



C(0-(-l)'«n- 



(-1)' 



(2.8) 
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with r| = 2. 

In contrast, the typical correlation function Ctyp (l) behaves 
quite differently. Its long-distance behavior is obtained by the 
following argument. Suppose spins S2 and ^3 are those to be 
decimated at a given SDRG step, as in Fig. 1. In that case, the 
correlation between S2 and equals — 3/4 + O (73/ £2)', while 
the correlation between ^4 and ^3 is of order — ^3 /£2. Thus, 
the typical value of the correlation function will be propor- 
tional to the typical value of // £2. Using the fixed-point dis- 
tribution (2.3), one finds ln|Ctyp(Z)| ^ — Ia//, i.e., the typical 
correlation function decays as a stretched exponential, with a 
nonuniversal prefactor i of order unity. 



B. Mean correlation function 

We now derive the mean correlation function in a more for- 
mal calculation which allows us to compute its amplitude in 
addition to its power-law decay. In the SDRG framework, the 
mean correlation function C{1) between spins separated by a 
distance / is obtained from the corresponding distribution of 
singlet-pair bond lengths in the ground state Ps{l), 



c(0 = --n(0, 



(2.9) 



since each singlet contributes a factor of —3/4 to C{1) and 
there are two spins in each singlet. 

The singlet-bond length distribution Ps{l) can be calculated 
from 

p,{i) = 2 «nP(y = n,/;n)dn , (2.10) 

Jo 

where P{J,i,Q.)dJdl is the probability of finding a coupling 
constant between J and J + dJ connecting spins separated by 
a distance between / and / + dZ at the energy scale £2, and the 



factor 2 comes from normalization. If we follow exactly the 
joint probability P (7, /; £2) along the SDRG flow, then we can 
obtain an exact expression for Ps{l)- In fact, we only need 
P(y,/;£2) at7 = £2. 

It turns out that we can carry out this task for A, = and 
couplings taken from the family of initial distributions. 



Po{J) = d{j)d{ao-j) 



1^0 /£2o 



£2() V / 



I --do 



(2.11) 



in which "do > gauges the strength of the initial disorder 
and £2o sets the initial energy scale. We first calculate the 
density of active spins «n- For that, we need P/(/;£2) = 
/ P(y,/;£2)dZ, which is obtained from the flow equation'** 



-^=P/(n;£2) J dJidJ3Pj{Ji;D.)PjiJ3;Q.) 

Introducing the Ansatz^^'^^ 

i}(£2) /£2N 



(2.12) 



PjiJ;^)^ 

into Eq. (2.12) yields 

i^(£2) 



£2 



1^0 



(2.13) 



(2.14) 



where F = In (£2o/£2). Thus, from the rate equation (2.4), we 
obtain 



1 



m = 



(l+i^oF) 



2 ' 



(2.15) 



and Eq. (2.5) is recovered in the low-energy limit F ^ oo. 

We now need to follow the SDRG flow of the joint distri- 
bution P (y, Z; £2), which is governed by the equation'^ 

dp r 

— d/id/2dZ3dyidy3P(/l,/i)P(£2,/2)^(^3,'3) 



5(/-/i-/2-/3)8 / 



£2 / 



(2.16) 



As shown in Appendix A, this can be done exactly by 
Laplace transforming P{J,l;Q.) and using an Ansatz for the 
corresponding flow equation. The final result for P{Q.,l) = 
P(y = £2,/;£2) is 



471 



n=l 



-(fro,(2.n, 



where a = i3-ov2'o, and Zq = 1 is the "bare" lattice spacing. 
Although the leading term of Eq. (2.17) had been obtained 
before,'^ the explicit dependence on the initial disorder distri- 
bution encoded in a was not emphasized. As will be shown 
next, this dependence is essential for our discussion. 
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Plugging Eqs. (2.15) and (2.17) into Eq. (2.10), we obtain 



87f 



n=l 



,-(f)-' 



(i+i^or)^r3 



-dr 



n=l 

3^,{l + o(VW^)}, 



\n+\ 



(-1)" 



2/o 



where 



1 

h Jo 



e ^6s. 



(V2/(7tn) + v/77^ 



(2.18) 
(2.19) 

(2.20) 



and we used/(Z > /o,n) ^ 1// |l + (^^/hT^ | in the last 
step. As explicitly shown in Eq. (2.18), the distribution of 
singlets in the ground state is independent of the initial disor- 
der parameter "do at all length scales. Moreover, it follows a 
universal power law in the large-distance limit. 

Finally, taking into account that singlets can only be formed 
between spins separated by distances corresponding to odd 
multiples of Iq, the mean correlation function takes the uni- 
versal form 



1 , if //Zq is odd, 
0, otherwise, 



(2.21) 



where X) = 1/4, irrespective of the initial disorder parameter. 
Note that Eq. (2.21) recovers Fisher's scaling result Eq. (2.8). 
In view of the fact that correlations between the spins in a 
singlet state are isotropic, correlations between components of 
the spins along a given direction a = x, y, or z should behave 
as 



1 



1, if ///o is odd, 
0, otherwise, 



(2.22) 



with a pref actor given by — X)/3 = —1/12. 

In order to check the prediction of Eq. (2.21), we calcu- 
lated the mean correlation function C{1) from numerical im- 
plementations of the SDRG algorithm on very large chains 
(2 X 10^ sites), with initial couplings following probability 
distributions of the form 



9(./-./mm)9(nQ-./) T^O /Ho 

1 - {Um/^of' 



(2.23) 



where Ho = 1, iJo > 0, and J^i^ > 0. Figure 2 shows, 
for various chains, the relative difference between the calcu- 
lated correlation function and the universal prediction, 5(/) = 
C{1) /Cu (/) — 1, as a function of /. We considered both the 
XX (A, = 0) and the isotropic Heisenberg (A,- =1) models, 
for which we averaged over 100 and 1000 disorder realiza- 
tions, respectively. In agreement with the previous analysis, 
the long-distance behavior of the correlation functions is well 
described by the universal prediction Cu{l), regardless of the 
model under consideration, within an error of less than 5%. 




Figure 2: (Color online) Relative difference 6 = C/Cu — 1 between 
the mean correlation function C and the universal prediction Cu as 
a function of the distance / in the SDRG framework, for various 
choices of initial disorder, and both XX and isotropic Heisenberg 
(or XXX) models. Chains A, D, F, and H have couplings distributed 
according to Pq {J) [see Eq. (2.23)] with (J^aim'^o) equal to (0.5, 1), 
(0,3), (0, 1), and (0,0.3), respectively. The results for chains Fxx 
and Hxx (omitted for clarity) are statistically indistinguishable from 
those for chain Dxx- Error bars (not shown for clarity) are of the 
order of the statistical data fluctuations. Lines are guides to the eyes. 



Moreover, the mean correlation function of chains Dxx, Fxx, 
and Hxx (all of which have Jmin = 0, as described in the fig- 
ure caption) are statistically identical at all length scales, in 
agreement with Eq. (2.18), which predicts the same short- 
distance behavior for those spin chains whose coupling con- 
stants are distributed according to Eq. (2.11). Notice that 
5(/) approaches zero for large I even for distributions with 
■^min > 0, which clearly do not belong to the particular class of 
distributions [Eq. (2. 1 1)] employed in the derivation of Cu (/). 

The clear difference between the convergence rates of the 
mean correlation functions in the XX and Heisenberg mod- 
els is due to the extra numerical prefactor of 1/2 present 
in the recursion relation of the latter model [cf. Eq. (2.1)], 
which delays the convergence of C(/) to the asymptotic form 
Cu (/ ) . This prefactor (which becomes negligible as the SDRG 
scheme proceeds) alters the relation between length and en- 
ergy scales, relevant for the derivation of Cu. However, at log- 
arithmically large energy scales, F = ln(i2o/i2) ^ In 2, the 
simple relation between length and energy scales in Eq. (2.7) 
is recovered. 



m. NUMERICAL RESULTS 

We now confront the predicted long-distance form of the 
mean correlation function, given in Sec. II, with numerical 
results for XX chains, obtained through the mapping to free 
fermions, and for isotropic Heisenberg chains, obtained by 
quantum Monte Carlo (QMC) calculations. 
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Figure 3; (Color online) Dependence of —1-0' (/) on the spin sepa- 
ration / in the XX model for three different probability distributions 
of the couplings: a box distribution (■&{) = 1) with J^[„ = 0, a box dis- 
tribution with y^in = 1 /4, and a binary distribution with y^j,, = 1/10 
[see Eqs. (3.1) and (3.2)]. The orange solid line corresponds to the 
disorder-free prediction (Ref. 4) —f'Cl'- = I/tt^. For short-length 
scales, all chains approach the behavior of the uniform system. After 
a disorder-dependent crossover length, the data approach the univer- 
sal prediction of Eq. (2.22), which is indicated by the black dashed 
line. Statistical fluctuations increase with /, and so results for / > 300, 
as well as error bars, are omitted for clarity. 



A. XX chains 

We analyzed disordered XX chains with periodic boundary 
conditions, and coupling constants following box distributions 



)Q{no-J) 1^0 

J 



(3.1) 



1 - (yn,i„/£2o)*» iio 
with fio = 1 , tS^o > 0, and Jmin > 0, or binary distributions 

Go (/) = ^8 (7 - 7„,i„) + is (7 - £2o) • (3.2) 

Below, we present results for different choices of parame- 
ters. Figure 3 shows the mean longitudinal correlation func- 
tion C"'{1) = {S'l'S^'^i^^ as a function of the spin separation / 
for a chain with 4000 sites and couplings taken from three 
probability distributions: two boxlike distributions (i3-o = 1) 
with /jnin = 1/4 and /,T,in = 0, and one binary distribution 
with /min = 1/10, in which we average over 700, 1000, and 
800 disorder realizations, respectively. Other disorder dis- 
tributions give similar results. The short-length behavior ap- 
proaches the uniform-system result,^ C" (Z) = — (nl)^^. After 
a disorder-dependent crossover length, I'^'^o the mean longitu- 
dinal correlation function decays as a power law with expo- 
nent Tj = 2, and the prefactor clearly approaches the universal 
value — 1/12 [see Eq. (2.22)]. Although not shown in the fig- 
ure, the typical longitudinal correlation function Cfj^p(Z), in 
contrast, has a nonuniversal prefactor 



Figure 4: (Color online) Dependence of |C" (/)| on the spin sep- 
aration / for the XX model with the same coupling distributions as 
in Fig. 3. The short-distance behavior, as in the uniform system, 
corresponds to a power law with exponent Tjc = 1/2. The random- 
singlet nature of the ground state governs the long-distance behavior, 
with r| = 2, but nonuniversal prefactors. The orange solid line cor- 
responds to the clean-system transverse correlation function (Ref. 5) 
(—1) 0.14709/\/7. The black dashed line is the prediction for the 
disordered system given in Eq. (2.22) (see text). 



It is remarkable how the random-singlet hallmark appears 
in Fig. 3. Since C^^ ~ l^^ both in the clean and in the dis- 
ordered cases, one could naively think that disorder does not 
play any role for C^^. However, our statistics are good enough 
to distinguish the different prefactors. We stress that the fluc- 
tuations seen at larger length scales reflect only the increas- 
ingly and inevitably poorer statistics, since the number of sin- 
glet pairs decreases as l^^ and their relative fraction becomes 
smaller and smaller Indeed, the mean correlation function is 
self-averaging. This could be directly double checked from 
the decrease of the relative fluctuations with the square root of 
the inverse chain size. 

We now turn our attention to the transverse mean correla- 
tion function C"(Z). Compared to its longitudinal counter- 
part, it behaves quite differently. The short-distance behav- 
ior, as in the uniform system, corresponds to a power-law de- 
cay with exponent Tjc = 1/2. After a nonuniversal crossover 
length, the random-singlet behavior is recovered, with a uni- 
versal exponent r| = 2, but nonuniversal prefactors, as shown 
in Fig. 4 and previously pointed out in Refs. 18,19,20. With- 
out loss of generality, the large-distance scaling form of C'"(/) 
can be written as 



C"(/»l) 



vll-\ 



if / is odd, 
otherwise. 



(3.3) 



with suitably chosen functions D'^ and v^. When couplings 
are drawn from disorder distributions sufficiently close to the 
IRFP form of Eq. (2.3), and Vg are expected to approach 
the constant values X)^ = — 1 /4 and = 0. 

The results from different coupling distributions provide 
evidence that and xjI indeed approach constant values for 
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Figure 5: (Color online) The summed transverse correlation func- 
tion Csu^(/) = C-"(/) +C-""(/+ 1) as a function of the distance 
/ (for odd /) for XX chains and the same coupling distributions 
as in the previous figure. Although for sufficiently strong disor- 
der (circles and triangles) the curves approach a power law, corre- 
sponding to the random- singlet exponent t] = 2 and to a prefactor 
("Do+D'e) /3 ~ —1/12, the less disordered system (squares) shows 
a different prefactor. Again, results for / > 400 and error bars are 
omitted for clarity. 



arbitrary initial disorder, with — and assuming close 
(but certainly distinct) values. Additionally, it seems that 
the quantity + Vg approaches an asymptotic value close 
to —1/4 for sufficiently strong disorder. This can be seen in 
Fig. 5, where we plot (for I odd) the combination C^^^ (/) = 
-[C"(/) +C"(Z + 1)]. Notice that, for the box distribution 
with y,j,in = and the binary distribution with Jy^m = 1/10, 
the curves for Qyni(/) are reasonably well described by the 
scaling form 1/ (12/^) in the long-distance limit. However, 
this is not the case for chains with couplings drawn from the 
box distribution with J^nm — 1 /4, at least up to the sizes stud- 
ied (/ = 1 000, not shown). Indeed, we argue in Sec. IV that 
deviations from that scaling form should be expected for the 
transverse correlations in XX chains. 

Finally, we report that we have considered also smaller 
chains (1000 sites) but with more disordered distributions 
(■Anin = 0, with iJo = 0.3 or i3-o = 0.6). For the sake of clarity, 
we have omitted their data in Figs. 3-5. The mean longitudi- 
nal correlation function is remarkably well described by the 
naive SDRG prediction (2.22). The mean transverse correla- 
tion function, on the other hand, is well described by Eq. (3.3) 
withuj + v^ w -1/4- 



B. XXX chains 

We now present QMC results obtained for the SU(2) sym- 
metric model, 

H^Y^J,SrSi+i , (3.4) 




Figure 6: (Color online) Correlation function (— l)'C^-'(/) in the 
ground state of isotropic random AF Heisenberg spin- 1/2 chains 
[Eq. (3.4)] of length Lo = 200 with disorder strength W = 0.75. The 
quantum Monte Carlo results were obtained at T// = 1.5x 10^^ and 
averaged over A'samples = 500 realizations. When the distance / be- 
tween spins is even (circles), the asymptotic regime is described by 
C'~{1) ~ 0.9/ (black, solid line) whereas for odd / (squares), the 
best fit gives C''(0 ~ -0.98//^ (red, dashed line). 

with the random AF couplings y,'s distributed according to the 
box distributions 

P{J)^^Q{j-Jil-W))Q{Jil+W)-j) . (3.5) 

The QMC algorithm we use is based on a stochastic series 
expansion of the partition function. ^^'^^ This is a finite tem- 
perature T technique which, in principle, allows access to 
ground-state properties, provided T is chosen to be much 
smaller than the finite-size gap of the system D. oc Lq~. As 
already discussed in several works (see, for instance, Refs. 20 
and 36,37,38), the ground-state properties in random spin sys- 
tems can be very hard to access because extremely small en- 
ergy correlations might develop between distant spins or spin 
clusters. For random finite chains, the dynamical exponent z 
is formally infinite since we expect exponentially small cou- 
plings to develop at large distances between spins, so that 
Q. o= exp(— -y/Lo). In order to accelerate the convergence to- 
ward the ground state, we used the p-doubling scheme^* and 
thus performed the QMC measurements at temperatures as 
small as 4 X 10^^ in units of J. We show in Figs. 6 and 7 
QMC results for the average spin-spin correlation function in 
the ground state, 

1 ^samples ^ Lq / 2 

^*samples ry=l ^ /=l 

where we perform disorder averaging over A^sampies indepen- 
dent random samples, as well as space averaging along the 
periodic chains. Note that the SU(2) symmetry of the Hamil- 
tonian ensures that C~^{1) = C^>'(/) = C"(/). 

As studied in great detail in Refs. 19 and 20, there is a 
crossover phenomenon which is governed by the localization 
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Figure 7: (Color online) Summed correlation function Cl^^iJ) = 
\C~~{1) +C~*'(Z+ 1)1 in the ground state of isotropic random AF 
Heisenberg spin- 1/2 chains [Eq. (3.4)] of length Lq = 200 with dis- 
order strength W = 0.75 (red squares) and Lq = 100 with VK = 1 
(orange diamonds). The quantum Monte Carlo results were obtained 
at r/7 = 1 .5 X lO"-'' for W = 0.75 and r/7 = 4 x 10"^ for IV = 1 and 
averaged over A'samples = 500 realizations for each disorder strength. 
The dashed line is the 1/ ( 12/~) prediction. 



length ^ of the corresponding one-dimensionalJordan-Wigner 
fermions with random hoppings. In order to be in the asymp- 
totic regime, i.e., in the vicinity of the IRFP, we have to look 
at system sizes Lo ^ ^. For the SU(2) symmetric case, ^ has 
been estimated to be ~ 20 for W = 0.75 and ~ 10 for W = l.^" 
Thus, in order to study the IRFP asymptotic regime, we study 
two different systems: W — 0.75 with Lq = 200 sites (see 
Fig. 6) and 'W = \ with Lq = 100 (see Fig. 7). In Fig. 6, we 
first clearly see the crossover behavior for distances / < 20 
and then the IRFP prediction, Eq. (2.8), recovered for larger 
separations. On the other hand, the prediction of Eq. (2.22) is 
not verified, and we confirm the observation already made for 
the XX case. Again, we can write 

^aa /, ^ 1 / Do^"^ if I is odd, 

C (/»i;)_-|^^^_2 otherwise, ^^'^^ 

where i)o and De are disorder-dependent prefactors. 

Nevertheless, the universality is recovered when looking 
at the sum of the prefactors (see Fig. 6) + De — —0.08 ~ 
— 1/12. As shown in Fig. 7, the quantity 

Cs'um(0 = -p-^(0+C-^-^a + l)] (3.8) 

seems to behave as 1 / (l2/2) for W = 0.75 and W = 1. 



C. Discussion 

The origin of the apparent universality of + De = — 1/4 
is not obvious. It is clear that the breakdown of the SDRG 
prediction (Uq = — 1 /4, = 0) must be related to the fact 
that a collection of singlet pairs is not an exact eigenstate of 



the Hamiltonian for any finite disorder Although the SDRG 
method becomes asymptotically exact at low energies, dec- 
imations involving spins connected by couplings of the or- 
der of the initial energy scale i2o unavoidably lead to signifi- 
cant eiTors due to the fact that the calculation is perturbative. 
Thus, instead of singlet pairs, these steps should really involve 
blocks of three or more neighboring spins, so that correla- 
tions "spread" over a few sites (whose number decreases as 
the strength of the initial disorder increases), forming clusters 
of correlated spins. 

Figure 8 shows the correlation functions = ( SqS'^j ) and 



Cf = 



SgSj^ between a reference spin So and the y-th neigh- 
boring spin to the right Sj as a function of j for XX chains 
with couplings drawn from boxlike distributions [panels (a) 
and (b)] and a power-law distribution [panel (c)]. For the par- 
ticular realizations in the figure, we find that, decimating the 
chains according to the SDRG scheme. So should couple to 
Sj* in a singlet pair, where j* = 97, 439, and 297 for panels 
(a), (b), and (c), respectively. Indeed, there is a pronounced 
peak at j*, as expected from SDRG. In addition. So also de- 
velops strong correlations with a few spins adjacent to Sj* and 
So. As expected for a localized phase, these contributions van- 
ish exponentially at larger distances. 

We now can define two spin clusters: The first one is com- 
posed by So and its neighbors such that the magnitude of the 
transverse correlation function between So and a spin belong- 
ing to that cluster is bigger than a certain cutoff, say, 10-^ the 
second cluster is analogous to the first one, but replacing So by 
S;* . Interestingly, the sum of all the longitudinal correlations 
between the spins in the first cluster and the spins in the sec- 
ond cluster approaches — 1/4. We thus say that the correlation 
between So and S^* "spreads" among the spins in the clusters, 
and each of them acts collectively as a single spin, leading to 
the universal result that Do + De = — 1 /4 at large length scales. 
However, such feature is not verified when we consider the 
transverse correlation function. As we show in the next sec- 
tion, this is related to the lack of total spin conservation in the 
transverse direction. 



IV. EXACTLY SOLVABLE MODEL 

The discussion at the end of Sec. Ill suggests that the forma- 
tion of spin clusters is responsible for the failure of the predic- 
tion of universal correlation functions, Eq. (2.21). According 
to the Marshall-Lieb-Mattis theorem,^^'***^ the ground state of 
a spin cluster with antiferromagnetic couplings is a singlet, if 
the number of spins is even, or a doublet, if the number of 
spins is odd. Thus, at low enough temperature, a cluster with 
an even number of spins does not contribute to the magnetic 
properties of the chain, while a cluster with an odd number of 
spins can be represented by an effective spin-1 /2 object. 

To gain insight into the origin of the apparent universality of 
Do + De observed in the numerical calculations of Sec. Ill, we 
now consider a chain that, at low energies, can be interpreted 
as being composed of a certain fraction e of "effective" spins 
and a fraction 1 — e of remaining "original" spins. Each effec- 
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Figure 8; (Color online) The magnitude of the transverse C*"^ and longitudinal C'~ correlation functions between So and Sj as a function of 
j (for j ^ 0). Main panels (a), (b), and (c) correspond to different samples drawn from different disorder distributions [see Eq. (3.1)] whose 
parameters (ymim'^^o) are (0.25, 1), (0, 1), and (0,0.7), respectively. Insets (i) and (ii) show C" and C'' between spins S/ and Sj as a function 
of / for different values of / between —3 and 3 and for values of / around j* . In the SDRG picture for these particular realizations. So and Sj» 
should couple in a singlet state where /* equals 97, 439, and 297 in panels (a), (b), and (c), respectively. The lines in the insets are guides to 
the eyes. 



live spin represents a cluster with an odd number of original 
spins. For simplicity, we assume here that each effective spin 
represents a cluster with only three original spins. In addi- 
tion, we locate the effective spin at the position corresponding 
to the middle spin of the underlying three-spin cluster. With 
this restriction, if the effective chain contains N spins (origi- 
nal and effective ones), there are, in fact, = .^V (1 + 2e) un- 
derlying original spins involved. Moreover, if there are ^— 1 
spins between a given spin pair in the effective chain, we say 
that I is the "effective" distance between the spins in that pair. 
[This is not the same as the renormalized distance defined in 
Eq. (2.2).] 

Now, we choose the couplings in the effective chain from a 
probability distribution like that in Eq. (2. 1 1), with tS-q <C 1, so 
that it is sufficiently close to the infinite-disorder fixed-point 
distribution and the occuiTence of "bad" decimations is highly 
improbable. With this choice of couplings, it follows from 
Eq. (2.19) that, in terms of the effective lengths I, the effective 
singlet distribution is given by 




(4.1) 



for odd I, while P, (/) = for even I. Therefore, the average 
number of singlets of effective length I in the effective chain 
is given by 

ff, (/) = (0 , (4.2) 

with / restricted to odd values. 

In order to obtain the ground-state correlations in the un- 
derlying chain, we have to determine the distribution of sin- 
glet lengths in terms of the underlying original distances I, and 
these will depend on how many effective spins are located be- 
tween the two spins in a given singlet. Let us consider a sin- 
glet formed between spins separated by an effective distance 
I. Hence, there are / — 1 intermediate spins, m of which we 
assume are effective ones. Note that there are three possible 
types of singlets (see Fig. 9): a pair of original spins (type 0), 



m-^^^\. . type 

• [m m m ] type 1 

(■ ■ "^ ^^ • • ~~ [" ■ ■ ] type 2 

// 

/ odd /o=l 



Figure 9: (Color online) The various types of singlets. Black circles 
denote isolated original spins. Blue squares denote original spins 
belonging to effective three-spin clusters. 

one original spin and one effective spin (type 1), and a pair 
of effective spins (type 2). If the singlet is of type 0, then 
the underlying distance / is given by Z = / + 2m; if the singlet 
is of type 1, then 1^1 + 2m + 1; and for singlets of type 2, 
1 = 1 + 2m + 2. We immediately conclude that, while singlets 
of types and 2 are associated with odd underlying distances 
I (since I is odd), singlets of type 1 are associated with even 
underlying distances. This leads to the appearance of corre- 
lations between spins separated by even distances, as in our 
numerical calculations, and in contrast to the assumption of 
the SDRG approach. 

For a singlet of length /, the number of intermediate ef- 
fective spins varies between and / — 1 . The probability of 
finding exactly m intermediate effective spins is given by 

(';^)b'"(i-s/-'-'". 

Thus, the numbers of singlets of types 0, 1, and 2 with under- 
lying length / are given, respectively, by 

A^s,o(/)-^^x(l-e)2/?(0 , (4.3) 
A^s,i(/) = ^^x2£(l-e)/?(/-l) , (4.4) 
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and 



A^s,2(/) = ^^xe2/?(;-2) , 



with 



^m_2;^ /'^-2m-l\ £"-(l-e/'""'' 
3,„t'oV ^ J {l-2mf 



(4.5) 



(4.6) 



In the limit of large /, the sum in R{1) can be written as an 
integral, which can be calculated by Laplace's method, using 
Stirling's approximation. The final result is 



R{l)^l[l+2t)l-^ + 0{l-') 



(4.7) 



Consequently, 



N,,Q{l)^^{\-zfr^ + o{r^) (for odd/), (4.8) 



A^u(0 = f^e(i-e)r2 + o(r3) 



and 



Ns,2 (l) 



N 



2,-2 



1 



-3\ 



(for even/), (4.9) 



(for odd/). (4.10) 



In order to calculate the ground-state correlations C"" (/) = 



(SfSfj^i), with a = X, y, z, let us focus on a fixed odd value of 
/. There are contributions to C"" (/) and C"" (/ + 1) coming 
from singlets of type 0, with underlying length / ; from singlets 
of type 1 , with lengths / — 1 and / + 1 ; and from singlets of 
type 2, with lengths / — 2, /, and 1 + 2. If we denote by c'^^^ 
(•-mid) "weight" of a spin in either end (in the middle) of a 
three-spin cluster to the a component of an effective spin (see 
Appendix B), we can combine all contributions to write (see 
Fig. 9) 



C"" (0 = - ^ {a^s,o (/) + [N^A (I 1) +A^s,i (/ + 1)] c; 
+ [N,,2{l-2)+N,,2il + 2)] (c«d)2 

+ N.,2iI)[2ic^Mf + iCi,f]} , 



end 



(4.11) 



and 

C""(/ + l) 



4-N 



{A^s,i(/ + 1) 



'mid 



+ 2[A^,,2(/)+A^,2(/ + 2)]c«,c«y} , (4.12) 
so that, to leading order in /, we have 

C"" (/) = - ^'"^ {(1 - e) V4£ (1 - £) 



'-end 



and 

C""(/ + l) 



l/-2{2e(l-8) + 4e2c«,}4 



mid" 



(4.13) 



(4.14) 



Note that the above results are significantly different from 
the bare SDRG results of Sec. 11, most notably in that the 
average correlation is, in general, not zero for even /. Both 
C"" (/) and C"" (/ + 1) decay with the random-singlet expo- 
nent r| = 2, but with different prefactors and X)" , respec- 
tively. However, we have 



^ K + - - ^{1 - [1 - (2c^„d + ^m>d)]e}' 



(4.15) 



For the XXZ chain, irrespective of the initial anisotropy A, 
the z component of the total spin is a good quantum number, 
assuming the value S'^^i = in the (singlet) ground state. This 

means that the sum of the ground-state correlations (^S'-S^j^ 

between a given spin ; and all other spins j in the chain is equal 
to — 1/4. Since this is also true for an effective spin, it follows 
that 2Cgjjjj + cjj^y = 1 (as can be easily verified explicitly; see 
Appendix B), and we must have 



1); 



1 



(4.16) 



irrespective of the concentration £ of effective spins (and thus 
of the initial disorder). Furthermore, in the Heisenberg limit 
(A = 1), owing to the SU(2) symmetry, we also have 



1 



(4.17) 



This last result, however, is not valid for A < 1 . In particular. 



in the XX limit, for which 2cl 



end 



'mid 



: 0.9142, we obtain 



-i(l -0.0858e)^ 



(4.18) 



yielding a weak dependence on e. 

For the isotropic Heisenberg chain and for C"(/), the an- 
alytical results derived in this section are in agreement with 
the numerical results in Sec. Ill, strengthening the conjecture 
of a universal behavior for sum of prefactors of the longitu- 
dinal ground-state correlations in random XXZ chains. The 
presence of larger effective-spin clusters (as typically occurs 
for weaker disorder; see Fig. 8) should not change the con- 
clusions of this section concerning the universality of 
since the sum of all the weights of the original spins belong- 
ing to an effective cluster is identically 1 when computed with 
respect to the symmetry axis. 



V. ENTANGLEMENT ENTROPY AND ITS RELATION TO 
THE CORRELATION FUNCTION 

In this section, we discuss the relation between the entan- 
glement entropy 5'(/) and the ground-state mean correlation 
function C(/). 

The entanglement entropy between two complementary 
subsystems A and B is given by 



S{1) = -TipA In Pa = -Trpg In 



(5.1) 
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Figure 10: (Color online) Ground state of the infinite-disordered AF 
spin- 1/2 chain. The entanglement entropy between the subsystem 
inside the box of length / and the rest of the chain is equal to the 
number of singlets shared by them. In this case, 5(/) = 5 x so, with 
.so being the entanglement entropy of a singlet pair when one of the 
spins is traced out. 



where / is the length of one of the subsystems, 

Pa^TibP = Y,Wb\pWb) (5.2) 
{' 

is the reduced density matrix obtained by tracing out the de- 
grees of freedom of subsystem B in the ground-state density 
matrix p = |(|)) ((|)|, and 1 is a set of states spanning the 
degrees of freedom of B (with a similar definition for pg). 

In the SDRG framework, the ground state of random XXZ 
chains is a collection of independent singlet pairs, i.e., 

Lo/2 

|c^> = |0,-> , (5.3) 

i=l 

where |0,) denotes the ;-th singlet pair and Lq is the total num- 
ber of spins in the chain (see Fig. 10). 

As the entanglement entropy between two spins in a singlet 
state is sq = In 2,^' the total entanglement entropy due to a 
given choice of subsystems A and B is equal to times the 
number of singlet pairs in which one spin belongs to A and 
the other one to B (see Fig. 10). Using this fact, Refael and 
Moore^' calculated the mean number of times that each bond 
is decimated, which is equivalent to the mean number of sin- 
glet lines crossing a given boundary. They found that the mean 
value of the entanglement entropy grows as (yln /) /3, with 
Y = In 2 being a universal number This is reminiscent of the 
entanglement entropy in conformally invariant (clean) one- 
dimensional quantum systems, which increases as (c In /) /3, 
where c is the central charge, a signature of the universality 
class of conformally invariant systems. ^^'^^ In the clean criti- 
cal XXZ chain, c = 1 . 

We now rederive the mean entanglement entropy by relat- 
ing it to the distribution of singlet lengths [see Eq. (2. 19)] and 
thus to the SDRG mean correlation function [see Eq. (2.21)]. 
By definition, the mean value of the entanglement entropy 
S (/) between a subsystem of length / and the rest of the chain 
is the sum of the entropies of all subsystems of length I, di- 
vided by Lq. (In a Lo-site chain with periodic boundary condi- 
tions, there are Lo different subsystems with the same length.) 
The contribution of a given singlet of length 1^ depends on the 
relation between 1^ and /. If 1^ > I, there are 21 different sub- 
systems of length / whose boundaries are crossed by the sin- 
glet line [see Fig. 11(a)]. Likewise, 2Zs different subsystems 
have their boundaries crossed by a singlet of length U < I [see 




Figure 1 1 : (Color online) Schematic entropy counting, (a) When the 
singlet length (in this case Is = 5) is greater than the subsystem length 
(1a = 3), there are different subsystems whose right boundaries are 
crossed by the singlet, (b) Otherwise, when the subsystem length Ib 
is greater than 1^ (in this case, 1b = 4 and = 3), there are 4 different 
subsystems whose left boundaries are crossed by the singlet. 

Fig. 11(b)]. Thus, 

2sn { ' ^'^ 1 

^ {k=\ k^+l J 

where A^s (^s) is the number of singlets with length /s in the 
ground state and 5o = In 2 is the contribution of a singlet pair 
to the entanglement entropy. 

As shown in Sec. II, A^s {h) is simply related to the corre- 
lation function by A^s (/j) = — 4LoC(/s) /3. Thus, substituting 
Eq. (2.21) into Eq. (5.4), and noting that C{k)=Q for even k, 
we obtain, for / ^ 1 and Lo °°, 

8 f ' 1 

S{l)^--so\J^kC{ls) + l E C{h)\ (5.5) 

2 ( 1 /■i-T 1 1 1 \ / 

= 3^" (2/. ^'^+27t,|.-'^J+^ ^'-'^ 

= l\nl + b, (5.7) 

in which 7=^0 = In 2, while b and 1} are nonuniversal con- 
stants that depend on the short-distance details of C(/). In 
this way, we recover the result obtained by Refael and Moore. 
Yet another derivation of the above result is presented in Ap- 
pendix C. 

Note that Eq. (5.5) relates the mean value of the entangle- 
ment entropy to the mean correlation function. This relation 
is valid only in the context of infinite-randomness spin chains, 
where both quantities are dominated by rare spin singlets. In 
AF spin- 1/2 chains without disorder, for instance, such rela- 
tion is no longer valid, and the correct expression is far from 
simple (though an efficient valence-bond approach can be de- 
veloped to study block entanglement properties'*^). 

Contrary to the naive universal form [Eq. (2.22)] of the 
ground-state correlation function, which is found not to hold 
when confronted with exact diagonalization or QMC calcula- 
tions, the universal prediction of Eq. (5.7) is fully supported 
by numerical results (see Ref. 43) and, as shown in Fig. 12, 
does not depend on the initial disorder strength. In view of the 
relation between these two quantities, revealed by Eq. (5.5), 
the arising question is how these seemingly contradictory re- 
sults can be reconciled. 

We address this question by looking at the entanglement 
entropy of the exactly solvable model of Sec. IV. The ground 
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Figure 12: (Color online) Entanglement entropy of random XX 
chains with Lq = 1000 sites, and disorder of the form of Eq. (3.5) 
with, from top to bottom, VK = (no disorder, open circles), IV = 0.25 
(cyan circles), IV = 0.5 (orange triangles), W = 0.75 (red squares), 
and W = 1 (green diamonds). These are exact diagonalization re- 
sults averaged over lO"* < ^samples < 10^ realizations. The black 
dashed lines are fits of the form 5(/) = (1/3) x In 2 x In / + /^ with & = 
1.048, 0.925, 0.876, 0.849 for W = 0.25, 0.5, 0.75, 1. In the inset, 
the constant b is also plotted versus 5" 1 , where 6^ = (lny)2 - (inT) , 
and naively fitted to = 0.812 + 0.03456/8. 

state of the model can be viewed as a collection of singlets of 
three different types (see Fig. 9). So, although Eq. (5.5) can 
no longer be used, the entanglement entropy is still related to 
the distributions of singlet lengths, in analogy with Eq. (5.4). 
However, we must remember that an effective spin represents 
a three-spin cluster. It can be easily shown that the entangle- 
ment entropies between a three-spin cluster and a single spin, 
as well as between two three-spin clusters, are also given by 
50 = In 2. However, since we have to average over all different 
possible subsystems of a given size, we must take into account 
situations in which one of the spins in a three-spin cluster lies 
in a different subsystem than the other two. When averaged 
over all subsystems and singlet types, these "internal" con- 
tributions, being only boundary effects, lead to an additional 
constant term and thus do not affect the scaling properties of 
the entanglement entropy. Explicitly, we have 

= -i{5o(/)+5i (0+52(0} , (5.8) 

in which So {I), Si (/), and 52 (/) are the average entanglement 
entropies due to singlets of types 0, 1, and 2, respectively. 
From Eq. (5.4), taking into account the "internal" contribu- 
tions bj (e), we can immediately write Sj (l), for j = 0, 1, and 
2, as 

r / Lo/2 ^ 

Sjil) = bj (8) +2^0 <^ L ^^N^-J ('^) + ^ L ('^) ' 

l/s=l k=l+\ J 

(5.9) 

with N^,j (/) given by Eqs. (4.8)-(4. 10). Bearing in mind that 
A^s,o (0 ™d (0 ^re zero for even I, while A^s.i (0 is zero 



for odd /, we can use the fact that 

A^s,o(4)+A^s,i(/s-l)+A^s,2(/s) = y/r^ + o(Z-^) , (5.10) 

for odd /s, to conclude from Eq. (5.8) that 

5(/) = |ln/ + /7(e) , (5.11) 

again with y = .^o = In 2, and an £-dependent constant b (e), as 
in Eq. (5.7). 

Although this exactly solvable model yields a nonuniversal 
ground-state correlation function C(/), the entanglement en- 
tropy 5(0 does follow the universal form derived by Refael 
and Moore. ^' This last result and the numerical confirmation 
of the universality of S{1) (see Ref. 43 and Fig. 12) suggest 
that a description of the ground state of random XXZ chains 
in terms of a collection of independent singlets remains valid 
at sufficiently large distances, provided we use the notion of 
effective spins already discussed in the previous sections. 

Finally, we should mention that the notion of the nonuniver- 
sality of the correlation amplitude due to high-energy small- 
scale details was first considered and investigated by Fisher 
and Young. Later, Refael and Moore^' realized that such de- 
tails only contribute to the inaccuracies in the location of the 
low-energy effective spins, which lead only to a surface term 
contribution to the entanglement entropy, as we have formally 
shown here. 



VI. STRUCTURE FACTOR 

In this Section, we compute the static structure factor 

s"-iq)=j^ L (5X) 

Lo-l 

=27t £ e-''?'/'oC""(0, (6.1) 
/=o 

which is straightforwardly related to the mean spin-spin time- 
independent correlation function C"" (/) and is directly mea- 
sured in neutron scattering experiments. Indeed, neutron scat- 
tering experiments probe the dynamical structure factor 

5"(C0,^) = — y e--/(>-'-)/'o r dte''"(Sfit)Sf{0)) , 

(6.2) 

which reduces to 5" (q) in Eq. (6.1) after an integration over 

CO. 

There are three noteworthy properties of 5" (q): 

i«(^)=5«(-^) , (6.3) 

L ^"(?)- J^i-o, (6.4) 
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Figure 13: (Color online) The longitudinal structure factor in the XX 
model for various disordered chains of lengths up to Lq = 4000. As 
can be seen, they are practically indistinguishable. Inset (i) high- 
lights 5' for small q in which the curves become somewhat distin- 
guishable near q = 0.06:1. Moreover, they slightly deviate from the 
clean system prediction = q (dotted line) tending to a universal 
form S~ = K\q\ (dashed line), with K = 7t^/12 (see text). Inset (ii) 
shows S~ near q = Tl/2, where disorder is irrelevant. See Eqs. (3.1) 
and (3.2) for the definition of the disorder parameters (here ilo = 1 
and -do = 0). 

where the sum is over the first Brillouin zone, and 

^"(9 = 0) = ^((Cf) , (6.5) 

where S^^ is the total spin along the a direction and (•••) 
means its expectation value on the ground state. Hence, 
S~ (0) = for the XX model and S"''^' (0) = for the isotropic 
case. Note that, in the continuum limit, Eq. (6.4) leads to 

We now show our numerical results on the static structure 
factor for the disordered chain in the XX and XXX models. 
We anticipate that all three properties [Eqs. (6.3)-(6.5)] are 
obeyed by our numerical results. 

A. XX model 

Figure 1 3 shows the longitudinal structure factor in the XX 
model for the clean system [in which S^{q) ~ \q\ (see in- 
set)] and various disordered chains. Because the longitudinal 
mean correlation function depends very weakly on the disor- 
der (only through the crossover length) and its universal am- 
plitude in the disordered case is very close to the clean-system 
value (1/12 in comparison to l/7t^; see Sees. II B and III A), 
the longitudinal structure factor is, for all practical purposes, 
universal. 

In general, due to spin conservation and the fact that 
(/) = for even I, S' equals and K at the points q = Q 
and q ~ K, respectively. Moreover, it can be shown that 
S~ has inversion symmetry around the point q — k/2, i.e.. 




J I I I L 
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Figure 14: (Color online) The transverse structure factor in the XX 
model for various disordered chains of lengths up to Lq = 4, 000. In- 
set (i) highlights S~ for small q, where ^-'■'(O) ^ 0.194, 0.159, and 
0.128 for the box distributions with Jj^ia = 1/4, /min = 0, and the 
binary distribution, respectively. In all cases, 5"* (q) —5'' (0) ~ |^| . 
Inset (ii) shows the behavior of near ^ = n. It follows the charac- 
teristic divergence of the clean system prediction (dotted line) until a 
crossover vector q^ above which it saturates to a constant (see text). 



S~ {n/2 + k) + S"^' {n/2 - k) = n, for -n/2 <k< k/2, which 
implies 5' (n/2) = k/2. 

Let us now consider the effects of disorder more closely. 
Because the behavior of S'' for small q is dominated by the 
large-Z behavior of the longitudinal correlation function, it fol- 
lows that S~{q<^l) ^ K\q\, where K = kHI/3 = 71^/12 w 
0.822 is a universal constant [see inset (i) of Fig. 13]. In 
the same way, because C"{1) ~ for even I, the behav- 
ior near ^ = 7i is also dominated by disorder. In this case, 
7t — 5~ ~ 7t) k\k — q\. This is verified by the data of 
Fig. 13, but not shown for clarity. Finally, because the Fourier 
series [Eq. (6.1)] at q near k/2 selects small values of /, the 
behavior near ^ = 7t/2 is dominated by the clean-system pre- 
diction, 3' (q ~ 7t/2) \q\ as shown in inset (ii) of Fig. 13. 
Note that all these arguments are valid because C" (/) = for 
even I. 

We now discuss the behavior of 5^ (see Fig. 14), in which 
disorder plays a more prominent role. In the absence of disor- 
der, Cf- (/) « _f^/z2/f+i/(2/f) ^ (_i)'/r//i/(2^)^ where K ^ 1 
is the Luttinger liquid parameter [see Eq. (1.2)], and F and 
F « 0.14709 (Ref. 5) are positive constants. Since the first 
term is monotonic, it dominates the structure factor for ^ ^ 1 . 
Hence, < 1) - 5c = 0) \q\^^^- The second (stag- 
gered) term gives a subdominant contribution ^ q^. Near the 
AF peak, q = K; however, the second term dominates, yield- 
ing a divergent contribution, namely, S'^ {q = K — e)^ \e\^^^^, 
for |e| < 1. 

Quenched disorder dramatically changes the above sce- 
nario. Rewriting the transverse correlation function 
[see Eq. (3.3)] as C"(/»l) ~ K + W^)5,,odd/ (3/^) + 
(-l)'\)g/ (3/^) (where 5/,odd = 1, if / is odd, and 5/.odd = 0, 
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Figure 15: The structure factor as a function of q for different disor- 
der parameters [see Eq. (3.5)] in tlie isotropic Heisenberg model for 
chains of lengths up to Lq = 200. Inset (i) shows 5 " for ^ <t; 1 and 
compares it with the (field-theoretical) clean system prediction 5 " 
K\q\, with K = \/2 being the Luttinger liquid parameter, and the 
disordered system prediction 5" — > k|^|, with universal K = 71^/12. 
Inset (ii) highlights 5" near the AF peak q~%. The dashed line 
is the clean-system prediction, namely, 5" — A.ln(l — g/lt), with 
X = —Tt/2 (see text). 

otherwise), it becomes clear that S '' {q 1) is dominated by 
the first term: S"" {q <^ 1) - [q = 0) ^ -n^ (uj + uj) \q\ /3, 
where — (Do + 'U'e) = 1/4 h- nonuniversal contributions [see 
inset (i) of Fig. 14]. Note that this result strongly relies on 
the fact that is nonzero; otherwise, S'' {q I) would 

be dominated by the short length scale contributions to C", 
which follow the clean-system prediction. 

The characteristic AF divergence near q = n is suppressed 
by disorder, as shown in inset (ii) of Fig. 14. For q > qc = 
In/Lc, where Lc is the crossover length above which the cor- 
relation function is dominated by the disorder, saturates 
to a constant approximately equal to 27iF^,^j 1 / Finally, 
because Eq. (6.4) has to be satisfied, the decrease of the AF 
peak is accompanied by its broadening as disorder increases 
[see inset (ii) of Fig. 14]. 



B. XXX Model 

We now turn our attention to the isotropic Heisenberg 
model (see Fig. 15). 

Similarly to the XX model, the universal features of the 
correlation amplitude yield a universal structure factor for q <C 
1: 5" (^) k|^| , Va, since +T)e = 1/4 is universal. By 
coincidence, the clean system prediction also scales linearly,^^ 
S^{q 1 ) K\q\, however, with a different slope K ~ 1/2 
which is the Luttinger liquid parameter As shown in inset (i) 
of Fig. 15, {q) is linear, but apparently with no universal 
slope. However, although we have statistical fluctuations and 
finite-size effects, a close inspection of our data for q < 0.07 
shows that K |^| fits better than K \q\. 



The logarithmic divergence at the q~% point is suppressed 
by disorder. As shown in inset (ii) of Fig. 15, fol- 
lows the clean-system divergence up to a crossover vector 
qc = 2k /Lc, where Lc is the crossover length above which 
the correlation exponent follows the long-distance predic- 
tion of the disordered system. Beyond qc, 5" saturates at a 
non-universal constant proportional to 27t^^i^j y/ln I /I. We 
note that the clean-system prediction depicted by the dashed 
line, 5" — {T^/2)ln{l — q/n), is actually the prediction 
of the Haldane-Shastry model,^^ which is a good approx- 
imation to the Heisenberg model for q < 137i/14.^^''*^ For 
q> 137t/14, 5" diverges as [— ln(l — ^/tt)]^/^, consistent with 
Q"- (-l)'^/ln7/Zfor/> 1. 

C. Discussion 

Summarizing, 5" (q) is peaked at ^ = 7t, Va, in both mod- 
els, while it approaches zero at ^ = 0, reflecting the antiferro- 
magnetic quasi-long-range order. Near q = n, the low-energy 
behavior of the structure factor is dominated by the short- 
length scale behavior of the spin-spin correlation function, and 
thus its scaling is determined by the physics of the clean sys- 
tem. However, the true divergence is completely suppressed 
by disorder, and the peak becomes shorter and broader On 
the other hand, the structure factor vanishes universally ^ \q\ 
for q <^ 1 [for 5^ in the XX model, one may consider the 
quantity {q) — S'^ {q = 0)] as a consequence of two fea- 
tures: (i) the power-law scaling of the mean spin-spin cor- 
relation function C""(Z) ^ l^^, with universal exponent'^'^^ 
r| = 2 and (ii) the amplitudes X)o and being different in 
magnitude. Moreover, due to the universal features of these 
amplitudes, 5 " (g') — 5 " (g' = 0) vanishes as K\q\, with a uni- 
versal K = Jt^/12, if a is a symmetry axis, and a nonuniversal 
K, otherwise. 

We now briefly discuss a controversy that has appeared in 
the literature. The dynamical structure factor 5"(co,^) was 
theoretically studied in Refs. 49 and 50 within the SDRG 
framework (and thus, in the small-co limit) and experimentally 
studied in Refs. 51,52,53, mainly by measuring the local dy- 
namical structure factor S (co) [obtained when one integrates 
5(tt),<?) over^] for the compound BaCu2(Sio.5Geo.5)207. Pre- 
viously, it had been thought that this compound is a good ex- 
perimental realization of the random AF spin- 1/2 chain with 
quenched bond randomness, since both the experimentally 
determined static magnetic susceptibility and local dynami- 
cal structure factor were found to be in good agreement with 
the strong-disorder theoretical predictions.^' However, further 
and more precise measurements appeared to be in contradic- 
tion with the strong-disorder theoretical scaling of S (co).^^'^^ 
Interestingly, the magnetic susceptibility measurements re- 
mained in agreement with the theoretical prediction for the 
disordered system.^^ This led to a puzzle. Thermodynamic al 
quantities seem to be dominated by the physics of the dis- 
ordered system, whereas spin correlations seem to show the 
clean system physics. 

We tentatively argue that the low-energy behavior of the 
quantity C05 (co) investigated experimentally may not be dom- 
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inated by the physics of the disordered system, even if the 
system itself is governed by a strongly disordered fixed point. 
The quantity S (co) is obtained from an integration over all val- 
ues of q. Therefore, it is dominated by the antiferromagnetic 
peak q = n. Such large momentum reflects the shortest length 
scales of the time-dependent correlation function, whose be- 
havior is expected to be dominated by the physics of the 
disorder-free system, as in the case of the time-independent 
correlation function. Hence, the experimental determination 
of the full ^-dependent 5" {q, Co) would be highly desirable. 

Recently, it was shown^"* that the spin- 1/2 compound 
MgTiOBOs displays a remarkable random-singlet signature 
for the magnetic susceptibility in a wide and accessible range 
of temperature. It would be interesting to perform neutron 
scattering experiments on this compound in order to check 
the predictions shown here for the static structure factor. 



experiments may be able to probe the universal features of 
those amplitudes. 

We also rederived the average ground-state entanglement 
entropy S{1), relating it to the mean correlation function. In 
contrast to the nonuniversality of Vo + De when considering 
C{1) along a nonsymmetry axis, the universal form of S{1) 
first derived by Refael and Moore^' was shown to hold for 
the exactly solvable model with effective spins, in agreement 
with the numerical data presented in Sec. V and previously 
elsewhere. 

Note added. Recently, we became aware of Ref. 55 where 
the dynamical structure factor is also studied and the contro- 
versy mentioned at the end of Sec. VI is considered. 
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VII. CONCLUSIONS 

In this paper, we revisited the ground-state properties of 
random-bond antiferromagnetic quantum spin-1 /2 chains us- 
ing analytical and numerical tools. We focused on the ques- 
tion of the universality of the spin-spin correlation function 
C(Z) and of the average entanglement entropy S{1), as func- 
tions of the distance I, as well as the connection between them. 

By following exactly the SDRG flow of a family of 
coupling-constant distributions, we showed that the SDRG 
approach predicts a fully universal power-law scaling form of 
the pair correlation function C(/), in which both the prefactor 
and the decay exponent are disorder independent. However, 
the SDRG prediction is strictly valid only in the limit of infi- 
nite randomness. Exact diagonalization and quantum Monte 
Carlo calculations on finite chains showed that this purported 
universality does not hold, except for the correlations along 
the symmetry axis in the XX model. Moreover, these numer- 
ical results reveal different correlation amplitudes for spins 
separated by odd and even distances, Vo and De, respectively. 
Nevertheless, we showed numerical evidence that the combi- 
nation Do +X)e, at least for the XXX model and for correlations 
along the symmetry axis in the XX model, is indeed universal 
and agrees with the SDRG prediction if we consider that spin 
clusters themselves develop singletlike correlations. In other 
words, the correlations of random singlets spread among spins 
in the clusters. 

As the average number of spins in a cluster depends on 
the details of the coupling-constant distribution, correlation- 
function prefactors are nonuniversal. However, the conser- 
vation law for the total spin component along the symme- 
try axis guarantees that correlations "spread" over all spins 
in the effective-spin singlets, leading to the universality of 
Do + De = — 1 /4. This hypothesis was confirmed analyti- 
cally in an exactly solvable model, in which a number of 
three-spin clusters were introduced into an infinite-disorder 
random-bond chain. Interestingly, the fact that Do and are 
different in magnitude has important experimental relevance: 
the small-(7 behavior of the structure factor is dominated by 
disorder. We have argued that ^-resolved neutron scattering 
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Appendix A: CALCULATION OF P{a, I) 

Laplace transforming Eq. (2.16) with respect to the length 
variable / yields 



^ = -P (£2, X) I dJ,dJ,P (7i , X) P (73 A) 8 - ^ 

(Al) 

whereP(y,X) = / exp(-A,/)P(y,/)dZ. We now substitute the 

a(X Q.) ^-^^^•^'> 
P(/,X) = ^^%^(") (A2) 



(A3) 
(A4) 



£2 \J J 
into Eq. (Al) and find a pair of equations. 



5f " = 



with the boundary conditions a(k = Q,Q.) = P (A, = 0,£2) 
iJ(£2). Since 



(A5) 
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we find the solutions^^ 

Poc + c2tanh(cr) 



P 



a 



c + potanh (cF) 



c cosh (cF) + po sinh (cF) 



(A6) 



(A7) 



where c = c (A-) is a constant of the flow, defined by c" = p^ — 
a^. Moreover, c is a real number since P > a, which can 
be shown by considering an explicit calculation of the mean 
distance between the "active" spins. 



dll J dJP{J,l) 

-lim— / dJP(J,X) 
a{X,Q.) 



lim i 



1 - 



P(X,i2) 



(A8) 



As Z > 0, Eq. (A8) ensures that a < p. 
The boundary conditions lead to c(k 

c2 (X < 1 ) = a'-X-" + O (r+i ) . From 



• 0) < 1, and so 



I = 



"4 '-^ 

X^Q X\ p 



lim , „- 
X^o 2A.p2 



2 2 

2^r 



/o(l+iJor)' 



(A9) 



one finds T = 1, the last step coming from the definition 
/"= Iq/iiq,, where Iq is the initial mean distance between the 
sites. For simplicity, we consider that initially all spins are 
uniformly separated, i.e., = fo = /o is the lattice spacing. 
The constant a thus equals ■^q^/TIq and depends on the two 
parameters of the problem, the initial length parameter /o and 
the initial disorder parameter "do- This is a consequence of the 
fact that the magnitude of a coupling constant J shared by two 
spins a distance / apart is correlated with I, and thus the joint 
distribution P{J,l) cannot be written as Pj (J) P/ (l) . 
We can finally obtain P(i2,/) by Laplace inverting P 



in the appropriate scaling limit, i.e. X ^ 0, F 
flA,i/2r ^ 0(1). Thus, Eqs. (A7) and (A6) become 



a = flVXcosech (^aVXTj 
P = fl\/Xcoth (aVXr\ , 



respectively, yielding 



p{a,i) 



47l2 



-ly+Si^exp 



«7t 



)l 



and 



(AlO) 
(All) 



(A12) 



Appendix B: RENORMALIZATION OF A THREE-SPIN 
CLUSTER 



by an effective spin- 1/2 object. If the three-spin Hamiltonian 
is given by 



(Bl) 



with y > and < A < 1, then the ground state is doubly 
degenerate. Thus, we can define an effective spin So such that, 
in the doublet subspace, the original spins are represented by 
"weights" defined by^^ 



r>a „Olr>a r>a „C(r>Ol r>a „0i 



with a = X, y, z and 



^end — "-1 



A+ VA^TS 



(B2) 



(B3) 



A + VA2T8 



^end — ^1 — '-S 



x,y x,v 

'^mid = '^2 ~ 



(a + v^2T8) 



2 ' 



(B4) 



(B5) 



iA(^A + VA^ 

"mid = ^2 ^ Z ' s 2 

1 + i A + VAZ + S] 



(B6) 



Appendix C: ANOTHER DERIVATION OF THE 
ENTANGLEMENT ENTROPY IN THE STRONG-DISORDER 
RENORMALIZATION-GROUP FRAMEWORK 

Following Refael and Moore, 2' we exactly calculate the 
mean number of times a given bond is decimated, which cor- 
responds to the mean number of singlet links Ms crossing a 
given point in the chain at the energy scale Q.. Averaging over 
all the sites in the lattice, this is simply the sum of the lengths 
of all bonds decimated until the energy scale Q., divided by 
the chain length: 

Ms(n)=/""dn/ dlnnP{Q.,l)l . (CI) 
Ja Jo 



Using the results in Eqs. (2.15) and (2.17) we find that 

1 



Ms (£2) 




1 



(C2) 
(C3) 



A cluster of three spins Si, S2, and S3, connected by anti- 
ferromagnetic XXZ couplings, can be replaced at low energies 



where we have used the relation between length and energy 
scales in Eq. (A9). Again, the reader should be aware of an 
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extra 1/2 pref actor when integrating over /, due to the fact that 
singlet lengths are restricted to odd multiples of Iq. Finally, 
because any subsystem has two boundaries, 

S{l)=2soM,^^lnl + b , (C4) 

in which y = io = In 2 is a universal constant. Note that the 
constant b = —1/3 + O ('^'^^) presented here has no physical 



meaning. Although its value does not depend on the initial 
disorder strength in this derivation, deviations from such value 
are expected due to the spin clusters crossing the boundaries 
(see Sec. V). 
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